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Abstract 

Anomalous dimensions of Wilson operators with large Lorentz spin scale logarithmically with 
the spin. Recent multi-loop QCD calculations of twist-two anomalous dimensions revealed the ex- 
istence of interesting structure of the subleading corrections suppressed by powers of the Lorentz 
spin. We argue that this structure is a manifestation of the 'self-tuning' property of the multi- 
loop anomalous dimensions - in a conformal gauge theory, the anomalous dimension of Wilson 
operators is a function of their conformal spin which is modified in higher loops by an amount 
proportional to the anomalous dimension. Making use of the parity property of this function and 
incorporating the beta-function contribution, we demonstrate the existence of (infinite number 
of) relations between subleading corrections to the twist-two anomalous dimensions in QCD and 
its supersymmetric extensions. They imply that the subleading corrections to the anomalous 
dimensions suppressed by odd powers of the conformal spin can be expressed in terms of the 
lower-loops corrections suppressed by smaller even powers of the spin. We show that these rela- 
tions hold true in QCD to all loops in the large (3q limit. In the M = 4 SYM theory, we employ 
the AdS/CFT correspondence to argue that the same relations survive in the strong coupling 
regime for higher-twist scalar operators dual to a folded string rotating on the AdSsxS 1 . 



to the memory of Alexander Nikolaevich VasiVev 



1 Unite Mixte de Recherche du CNRS (UMR 8627). 



1. Introduction 



Success of QCD as the theory of strong interactions relies in part on the possibility to predict 
the scale dependence of various physical observables. In the most prominent example of deeply 
inelastic scattering, this dependence is governed by anomalous dimensions of composite Wilson 
operators that arise in the operator product expansion of conserved currents. The anomalous 
dimensions depend on the coupling constant as well as on the quantum numbers of operators 
(Lorentz spin, isotopic charge) and their properties reflect the symmetries of the underlying 
gauge theory. In the maximally supersymmetric M = 4 Yang-Mills theory, the gauge/string 
correspondence allows one to relate the anomalous dimensions at strong coupling to the energies 
of strings propagating on the AdSsxS 5 background [JJ. In the case of QCD, due to lack of 
a similar dual stringy description, one can only calculate the anomalous dimensions at weak 
coupling as series in the coupling constant. 

At present, the most advanced QCD calculations of anomalous dimensions have been per- 
formed to three-loop accuracy for the Wilson operators of twist two j2j [3]. The resulting ex- 
pressions for multi-loop anomalous dimensions are complicated functions of the Lorentz spin N 
carried by the Wilson operators but they simplify significantly for large N. In this limit, for 
N ^> 1, the anomalous dimensions scale logarithmically j(N) ~ IniV with subleading correc- 
tions running in inverse powers of N [4]. The explicit three-loop calculation [21 [3] revealed the 
existence of interesting structure of the subleading corrections. It turned out that, to three-loop 
accuracy, correction ~ 1/N to the twist-two anomalous dimension can be expressed in terms 
of the leading ~ iV correction evaluated to two loops only (MVV relation). Recently, it was 
shown in Refs. [5] that this relation can be explained by modifying the evolution (renormalization 
group) equation in such a way [6] as to preserve the reciprocity relation between its solutions 
(see Eq. flO} below). 

In this paper, we shall argue that the MVV relation is a manifestation of the following 'self- 
tuning' property of the multi-loop anomalous dimensions [7] - in a conformal gauge theory, the 
anomalous dimension is a function of the conformal spin of the Wilson operator which is modified 
at higher loops by an amount proportional to the anomalous dimension. Making use of the parity 
property of this function and incorporating the beta-function contribution, we shall demonstrate 
the existence of (infinite number of) MVV like relations between subleading corrections to the 
twist-two anomalous dimensions at large N in QCD and its supersymmetric extensions. We shall 
show that subleading corrections to the anomalous dimensions suppressed by odd powers of 1/N 
are not independent and can be expressed in terms of coefficients accompanying smaller even 
powers of 1/N. In the M = 4 SYM theory, we shall employ the AdS/CFT correspondence to 
argue that the same relations survive in the strong coupling regime for scalar operators of higher 
twist. 

The anomalous dimensions of twist-two operators govern the scale dependence of nonper- 
turbative parton (quark, gluon) distribution functions in the deeply inelastic scattering (DIS) 
[H [91 [10]. They are related to the parton space-like (s) splitting functions Ps(z) through the 
Mellin transform 



Jo 

The function Ps{z) defined in this way has a clear physical meaning. It defines the probability 
for a parton to split into another parton carrying the momentum fraction One encoun- 

1 The splitting functions cease to have a probabilistic interpretation at higher twists 




(1.1) 
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ters another set of nonperturbative parton fragmentation functions in hadron production in 
the e + e~ —annihilation which is just the cross-channel of the DIS process (see review [12] and 
references therein). The scale dependence of the fragmentation functions is governed by the 
time-like (t) splitting functions Pt(z) whose moments provide yet another anomalous dimen- 
sions 7T (iV) = - ^dzz N ~ 1 P T {z). In distinction with the parton distribution functions, the 
fragmentation functions do not admit the operator product expansion and, as a consequence, 
their anomalous dimensions Jt(N) are not related to composite local Wilson operators. The two 
functions are related however by the crossing symmetry. This suggests to write down the following 
reciprocity relations between the space-like, Ps{z), and time-like, Pt{z), splitting functions 

Fr(z)=P s {z), P T (z) = -zP s {l/z), (1.2) 

the so-called Gribov-Lipatov [13] and Drell-Levy-Yan [2] relations. Here in the first relation the 
two functions are equated within their physical region < z < 1, while in the second relation the 
space-like splitting function is analytically continued through z — 1. In general, the two relations 
in (11.2p are independent on each other. However, if they were fulfilled simultaneously, one would 
obtain a 'strong' reciprocity relation for P(z) = Pt{z) = Ps{z) 

P(z) = -zP(l/z). (1.3) 

It is known [T5J [T6], [TTJ, [18], [19] that this relation is satisfied to one loop only, while both relations 
in (11.21) are violated to two loops. 

Another set of relations for the twist-two anomalous dimensions comes from a space-time 
picture of parton (quark, gluon) splitting in the space-like (distribution) and time-like (frag- 
mentation) processes. In the latter case, imposing the condition of the angular ordering for 
successive partonic splittings one obtains the following remarkable functional equation for the 
leading asymptotic behaviour of the time-like gluon anomalous dimension 7t(-/V) at small N [21] 

a N 

1t(N) = ls (N - 7 t(A0) , ls(N) = + . . . , (1.4) 

where ellipses stand for subleading terms as iV — ■> and a s /N 2 = g 2 /(kirN 2 ) = fixed. The rela- 
tion (jl.4p leads to a quadratic equation for Jt(N) whose solution resums the leading corrections 
~ a™/N 2n ~ l to 7r(A r ) for all n > 1. It has been proposed in Refs. [21 E] that a relation similar 
to (11.41) should also hold at large N [f] 

-y*(N)=f(N-la 7<T (N)), (1.5) 

where a = — for the space-like and a = + for the time-like anomalous dimensions, 7_(iV) = 
75 (N) and j+(N) = 7^(A r ). The function f(N) is assumed to be the same for a = ± and, most 
importantly, its large iV expansion takes the form 

f(N) = A [if)(N + 1) + 7e] + B + • N- 1 + 0{N~ 2 ) , (1.6) 

2 Notice that the second relation in ()1.2|) still holds true to two loops for the so-called physical evolution kernels 
governing the scale dependence of the DIS structure functions and the hadron fragmentation functions [20] . 

3 Due to different conventions for the anomalous dimension (see Eq. (|2.ip below), "j a {N) differs from similar 
expression in Ref. [5] by the overall factor (—1/2). 
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where ip(x) = d\a.Y{x)/dx and 7 E is the Euler constant. Then, the relation ( 11. 5ft leads to the 
following expression for the twist-two (nonsinglet quark and diagonal elements of the singlet 
mixing matrix) anomalous dimensions at large A 

7(T (A) = AlnN + B + C a N- 1 \nN+ (A, + |A) A" 1 + e>(A" 2 ), (1.7) 

with TV = Ae 7E . Here the leading term is related to the universal cusp anomalous dimension [4] 

A = 2T cusp (a s ), (1.8) 

the B— coefficient is the same for the space-like and time-like distributions but it depends on 
the quantum numbers of partons (flavour, spin projection). From (11.51) and (11.61) one finds the 
remaining coefficients as [5] 

^w) = _i aA 2 {as) ; Da = -l aA {a s )B{a s ) . (1.9) 

The expression for C a is in agreement with explicit three-loop QCD result for the space-like [21 [3] 
and (nonsinglet) time-like anomalous dimensions [22], while for D a the mismatch is proportional 
to the beta-function 

pCMW) = -± a A(a s )B(a s ) - \A{a s )(3{a s ) . (1.10) 

As was emphasized in Refs. [21 E], this calls for a further structural explanation of subleading 
1/N corrections to the anomalous dimensions (I1.7P at large A and, more generally, for a better 
understanding of the origin of the functional relations (11.4|) and (11.51) . 

As we will see momentarily, for the space-like anomalous dimensions 7s(A), the relation 
( 11.51) is ultimately tied to conformal invariance of a classical four- dimensional Yang-Mills theory. 
Moreover, the conformal invariance allows one to extend the relation (ll.5p to anomalous dimen- 
sions of the so-called quasipartonic operators [23] of arbitrary twist L. A distinguished feature 
of these operators is that they are built from exactly L fundamental fields of twist 1 and from an 
arbitrary number of covariant derivatives projected onto the light-cone. The twist-two operators 
correspond to L = 2. 

The classical Yang-Mills Lagrangian is invariant under conformal transformations but this 
symmetry is broken on the quantum level unless the beta-function vanishes to all loops (see 
review [21] and references therein). If the conformal symmetry was exact (as in M = 4 SYM the- 
ory), the quasipartonic operators could be classified according to representations of the collinear 
SX(2;R) subgroup of the full 50(2,4) conformal group [22]. Conformal invariance ensures that 
the operators belonging to different SL(2; M) multiplets cannot mix under renormalization and 
their anomalous dimension, 75(A), depends on the conformal SL(2; R) spin of the multiplet. By 
definition, the conformal spin of the quasipartonic operator of twist L equals j = ^(N + £) with 
and i being its Lorentz spin and the scaling dimension, respectively. The former is protected 
from perturbative corrections, whereas the latter equals N + L to the lowest order in the coupling 
constant and receives anomalous contribution 75(A) at higher orders. As a result, the conformal 
spin of quasipartonic operators gets modified in higher loops as [261, 123 U\ 

Jba re = A + ±L ^ j = N+lL+±y s (N). (1.11) 

Then, the conformal symmetry implies that the anomalous dimension of the quasipartonic oper- 
ator is a function of j, or equivalently 

75(A) = (A + | 7 5(A)), (1.12) 
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where the scaling function f(N) depends on twist L and other quantum numbers of the operator. 
Here the superscript indicates that the function /(^ =0 ) is defined in gauge theory with vanishing 
beta-function. The relation (11.121) coincides with the first relation in (11.51) in the space-like case 
a — —. 

In gauge theory with nonvanishing beta-function, the relation (11.121) should be modified to 
incorporate the additional conformal symmetry breaking corrections. The latter are known to 
be renormalization scheme dependent. The analysis can be simplified by renormalizing the 
Wilson operators in the dimensional regularization scheme (DREG) with d = 4 — 2e% In this 
case, the coupling constant acquires a nonvanishing dimension and the beta-function in (4 — 
2s) — dimensional gauge theory is given by 

(3 £ {a s ) = -2e + (3{a s ) . (1.13) 

A crucial observation is that (3 £ {a s ) vanishes at e cv = fi(a s )/2 and, therefore, the gauge theory is 
conformal in d cr = 4 — 2e CT dimensions. This opens up the possibility to calculate the anomalous 
dimensions within the e— expansion using the conformal field theory technique as pioneered in 
Refs. [28], [29]. As before, the quasipartonic operators can be classified according to representation 
of the SL(2; R) collinear subgroup of the conformal group in d CT — dimensions and their anomalous 
dimension is a function of the conformal spin j cr = ^(N + £ CT ). The only difference as compared 
with the d — 4 case is that for e CT ^ the scaling dimensions of fundamental fields get shifted by 
(— e cr ) so that the scaling dimension of the quasipartonic operator built from L fields takes the 
form £ CT = N + L — Le cr + , ~fs(N). Then, the conformal invariance in d cr — dimensions implies that 
the anomalous dimension of twist— L operators in d = 4— dimensional gauge theory is a function 
of the conformal spin j CT 

ls(N) = (N + ± 7 s(A0 - \L(3{a s )) . (1.14) 

It is important to stress that, in distinction with f^ =0 \N), the function f^ Q \N) depends on 
the renormalization scheme. Invoking the same arguments that led to (]1.5j) . one can extend the 
relation (11.14j) to the time-like anomalous dimension 

1t{N) = fW°) (AT - \ lT {N) - \ma s )) . (1.15) 

The equations (j!.14p and f 1 1 . 1 5 j) generalize the relations (jl.5p proposed in 0, [5]. They relate 
the space-like and time-like anomalous dimension of twist L to the same function f^°\N) and 
incorporate the beta-function contribution. 

Conformal symmetry allows one to relate the space- like anomalous dimension 75 (N) to yet 
another function f(N) but it does not tell us much about the properties of the latter. Therefore, 
to make use of the functional relation (11.141) it has to be supplemented with additional relation 
for the scaling function f(N). As a hint, let us consider the twist-two anomalous dimension (11.71) 
and examine the first few terms in the large N expansion of f(N), Eq. (II. 6p . Their substitution 
into (11.141) and (I1.15P yields (for L = 2) the large N expansion of the space-like and time-like 
anomalous dimensions (11.7j) with the D— coefficient matching the exact three-loop result (ll.lOp . 
Notice that it is due to zero value of the coefficient in front of iV -1 in the right-hand side of (11.61) 
that the C— and D— coefficients parameterizing subleading ~ N^ 1 corrections to the anomalous 
dimension can be expressed in terms of the leading coefficients, Eqs. (11.91) and (11.101) . The 



SYM theories, to preserve the supersymmetry, one uses the dimensional reduction scheme (DRED) instead. 
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question arises whether similar relations also exist for the subleading ~ N~ n corrections (for 
n > 2). We shall argue that the answer is positive and it relies on the following remarkable 
property of the function f(N): the large N expansion of the function f(N) only runs in integer 
negative powers of the parameter 

J 2 = N(N + 1), (1.16) 

which has the meaning of a 'bare' quadratic Casimir of the collinear SL(2; R) group. For twist- 
two operators, the same property can be expressed as (with J = Je 7E ) 

f(N) = r cusp K)in J 2 + /(°) + f^/j 2 + ;( 2 )/j 4 + oil/ j 6 ) , (i.i7) 

where f^ k > (with k > 0) are given by series in a s with the coefficients depending on In J. The 
relation (11.17j) implies that each term of the asymptotic series (11. 17ft is invariant (modulo con- 
tribution from In J) under J — > —J, or equivalently N — > — 1 — N. In what follows we shall 
refer to (I1.17P as a parity preserving asymptotic series. We will show that the relations (I1.14[) . 
(I1.15P and (I1.17P are in a perfect agreement with available two- and three-loop results for various 
anomalous dimensions in QCD and in supersymmetric Yang-Mills theories. 

The relation ( 11. 17jl has a number of important consequences. Firstly, going back from the 
N— space to the z— representation, f(N) = — dz z N ~ 1 Pf(z), one finds that (11. 17ft is translated 
into the reciprocity condition (I1.3P for the splitting function Pf(z). Secondly, substituting (11.171) 
into (11.141) and (11.151) . one obtains that the coefficients in front of N~ 2n and N~ 2n ~ l in the 
large N expansion of the twist-two anomalous dimensions are related to the same coefficients 
with < m < n. Excluding the latter, one gets an infinite set of relations between the 
former (odd and even) coefficients, including the relations (II. 9p and (11.101) for n — 0. Moreover, 
re-expanding the anomalous dimensions at large TV in inverse powers of J = y/N(N + l), one 
finds that corrections suppressed by odd powers of J can be expressed in terms of the coefficients 
accompanying smaller even powers of J to less number of loops. 

The relations ( 11.141) and ( 11. 17ft are universal properties of the twist-two operators in gauge 
theories - they hold true both in QCD and in supersymmetric Yang-Mills theories with even 
number of supercharges M = 0, 2, 4, to three loops at least. In M = 1 theory the function f(N) 
develops an additional square-root singularity ~ a/1/ J 2 due to degeneracy of the parity-even and 
parity-odd anomalous dimensions in the large iV limit. Finally, the relation (11.1 7[) can be also 
tested to all loops: In QCD, we will determine the function f(N) in the so-called large (3q— limit 
and show that it verifies (I1.17P indeed. In the M = 4 theory, we will apply the gauge/string 
correspondence to evaluate f(N) in the strong coupling regime for the special class of scalar 
operators of higher twist dual to a folded string rotating on the AdSsxS 1 part of the target space 
of the type IIB string theory [301 EI]- 

The paper is organized as follows. In Sect. 2, we apply well-known one- loop expressions 
for various twist-two anomalous dimensions in QCD and in SYM theories to justify the parity 
preserving relation (I1.17P and to establish its connection with the reciprocity relation (ll.3p . In 
Sect. 3, we use explicit results for two- loop and three-loop anomalous dimensions available in the 
literature to verify the relations (11.141) . (11.151) and (I1.17P beyond the leading order. Then, we 
make use of the parity preserving property (I1.17P to find the relations between different terms 
in the large N expansion of the anomalous dimensions analogous to (II. 9p and (11.101) . In Sect. 4, 

5 The same property would not be manifest if one replaced J 2 by its expression (|1.16j) and re-expanded f(N) 
in powers of 1/N. The resulting series for f(N) involves both even and odd powers of 1/N and its first three 
terms match (|1.6p for /(°) = B. 
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we argue that the property (11.171) holds true in QCD for twist-two anomalous dimensions to all 
loops in the large (3q limit. In Sect. 5, we employ the AdS/CFT correspondence to determine 
the function f(N) in the planar M = 4 theory in the strong coupling regime for scalar operators 
of higher twist. Section 6 contains concluding remarks. Some details of our calculations are 
presented in the Appendix. 



2. Symmetries of twist- two anomalous dimensions 

Let us start with the definition of the space-like anomalous dimensions in a generic four-dimensi- 
onal Yang-Mills theory with the SU(N C ) gauge group. In general, the Wilson operators with the 
same Lorentz spin iV mix under renormalization and verify the Callan-Symanzik equation 

mJ-r)o { °\0) = -t*(N)O®(0), (2.1) 



d In fx d In A 

where the coupling constant A is defined slightly differently in QCD and in SYM theory 

\ (QCD) _ 9s _ ^_S_ \(SYM) _ 9~N C . . 

~ 8tt 2 ~ 2vr ' 8tt 2 ' 1 ' ' 

Here the mixing matrix is given by a series in the coupling constant 

7(A0 = A 7 o(A0 + A 2 7i(A0 + -- - , (2-3) 
with 7fc(iV) being matrices, and the beta- function is defined as 

^ = /3(A) = -/3oX - PiX 2 + 0(A 3 ) . (2.4) 

The beta-function coefficients are given in QCD by 

$ QCD) = y N c - 2 -n f , /^ CD) = y Nl - jj N c n } - C F n f (2.5) 

with rif the number of quark flavours and Cp = (N 2 — 1)/(2N C ) the quadratic Casimir in the 
fundamental representation of the SU (N c ) group. In the SYM theory with Af supercharges one 
has 

^(sym) = (4 _ jy) ? ^(svm) = {2 _ M) ^_ M y (2.6) 

For M = 2 the exact beta-function is given by the one-loop expression while for M = 4 it vanishes 
to all loops and the corresponding gauge theory is conformal. 



2.1. Nonsinglet anomalous dimensions in QCD 

Let us first consider nonsinglet twist-two quark operators in QCD. They take the form 

^ ns) (0) = q a>i T? 3 r^f-^O) , (2.7) 

where q^j denotes the quark field component carrying the flavour j and the spinor index (3. 
Also, D + = d + + igA + is the light-cone component of the covariant derivative. The matrices 
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T°j stand for the generators of the flavour SU(n/) group while the spinor matrices T a p serve to 
select the so-called "good" (twist-one) components of quark and antiquark fields either with a 
definite helicity, or with a transverse polarization [21]. A distinguished feature of the operators 
( 12.71) is that they can not mix under renormalization with other operators but only with the 
operators d™O^ n (0) containing total derivatives. This mixing can be avoided by going over 
from the operators (0) to their forward matrix elements with respect to some reference state 

<^ ns) (o)>. 



Depending on the choice of the T— matrices in (12. 7h one can distinguish three different twist- 
two operators: unpolarized (r = 7+), longitudinally polarized (r = 7+75) and transversity 
(T = 7+7^) operators|§ To the lowest order in the coupling constant their anomalous dimensions 
are given by well-known expressions [23] 



7o(A) = C f 



4^(iV + l)+4 7E 



1+r) 



N(N + 1) 



(2.1 



where 77 = 1 for the unpolarized and longitudinally polarized operators and rj = — 1 for the 
transversity operator. Here ip(x) = d\nT(x)/dx and 7 E = — is the Euler constant. In the 
momentum fraction representation (11.11) . the corresponding splitting functions read 



Po(z) = C F 



(1-*)+ 



+ 35(1 -z) + (l-z) (1+r/) 



where the '+' distribution is defined in a standard way 



(!-«)+ 



1 - z 



(2.9) 



(2.10) 



with ip(z) being an arbitrary test function. 

As was already mentioned in the Introduction, to one-loop order the space-like and time-like 
splitting functions satisfy the reciprocity relations (11.21) . Indeed, one verifies that for z ^ 1 the 
splitting function Pq(z), Eq. ( 12. 91) . fulfills the relation ( 11. 31) . The same property in the A— space 
finds its manifestation in the asymptotic behaviour of the anomalous dimension 70 (A) for large 
A. In this limit, the anomalous dimension (12.81) scales as 70 (A") = 4CVln A + . . . with subleading 
corrections running in powers of 1/A. However, changing the expansion parameter from A to 
J 2 = A(A + 1) and expanding the anomalous dimensions in inverse powers of 1/J one finds that 
the subleading corrections run in even powers of 1/ J only 



70(A) = C F [2 In P - 3 - (± + rj) J~ 2 - ± J~ 4 + ^ J" 6 + O (J^ 8 )] 



(2.1i; 



where J = Je 7E . As follows from (11.141) . the functions 75 (A) and /(A) coincide to the lowest 
order in the coupling constant, /(A) = A7 (A) + 0(X 2 ), but differ from each other starting 
from two loops. Matching (12. lip into (11.171) . one can identify one- loop contribution to the cusp 
anomalous dimension r cusp (a s ) and to the coefficients f( n \ 

Each term of the expansion (12.111) is invariant under substitution J — ► — J, or equivalently 
A — > — 1 — A, whereas the exact expression (12.81) is not 



7o(A) - 7o(-A - 1) = -4CV7rcot(7rA) . 



(2.12) 



6 In the case of transversity, the operators (j2.7p evolve autonomously even for flavour singlet quark distributions. 
The same holds true for the linearly polarized gluon distribution |32j . 
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Still, the two relations, Eqs. ( 12.111) and ( 12.121) . are consistent with each other - the one-loop 
anomalous dimension ( 12.81) has (Regge) poles at integer negative N and expansion ( 12.111) is only 
well-defined for Re N > 0. Regge singularities of ~fo(N) lead to a factorial growth of the expansion 
coefficients in the right-hand side of (12.111) (they scale as ~ (— 1 ) n (2n)\/(n(2-nJ) 2n ) for n > 1). 
Let us now show that the validity of the reciprocity relation (II. 3p for the splitting function and 
the absence of odd powers of 1/ J in the asymptotic expansion of the anomalous dimension ( 11. II) 
are in the one-to-one correspondence with each other. To this end, one changes the integration 

1 /2 

variable in ( 11.11) as z = e~ X//j with j = N + \ — (J 2 + |) and evaluates derivative 7'(iV) = 
<9yv7(iV) as 

1 f°° 

j'(N) = — dx e~ x [ie l/(2j) P(e- x/j )] . (2.13) 
j Jo 

Expanding the integrand at large j and integrating term-by-term one can obtain an asymptotic 
expansion of 'y'(N) in powers of If asymptotic expansion of 7(iV) run in even powers of 1/J, 
or equivalently 1/j, then similar expansion of the derivative y(N) would run in odd powers of 
1/j leading to 

xe */m p(e~ x/j ) = -xe~ x/{2j) P(e x/j ) . (2.14) 

For x > and z = e~ x ^ this relation is equivalent to the reciprocity condition (II .3p . In the 
similar manner, if the splitting function verified the reciprocity relation (II. 3p . the corresponding 
anomalous dimension would be given at large N by the parity preserving asymptotic series (j!.17p . 

Explicit calculation shows [TBI [TT] that for the nonsinglet twist-two quark operators the reci- 
procity relation (11.31) breaks down starting from two loops and, as a consequence, the asymptotic 
expansion of the nonsinglet anomalous dimension involves odd powers of 1/ J at higher loops. We 
will demonstrate in Sect. 3 that the function f(N) related to the nonsinglet anomalous dimension 
through (11.141) still possesses the parity preserving property (11.171) to three loops at least. 

2.2. Singlet anomalous dimensions in QCD 

As the next step, we consider flavour singlet quark and gluon operators of twist two 

O% ) (0) = q a>i r af) D»- 1 qp, i (0), O$(0) = tr D»~ 2 F v+ (0)] , (2.15) 

where = F® + t a = -{D^, D + ] is the light-like component of the gauge field strength tensor 
and t a are the generators of the SU(N C ) gauge group in the fundamental (quark) representation. 
As before, the spinor matrix T a p and the Lorentz tensor IL^ select 'good' (twist-one) components 
of the quark and gauge fields. For T = 7 + and IT^ = g^ u the matrix elements of the operators 
( I2.15P define unpolarized quark and gluon distributions, while for Y = 7+75 and 11^ = e M ^ + _ 
they define longitudinally polarized distributions. 

The operators (12.151) mix under renormalization with each other and with operators involving 
total derivatives. As before, the latter can be eliminated by going over to the forward matrix 
elements in ( I2.15p . Then, the matrix elements (0^(0)} and (0^(0)) verify the evolution equa- 
tion (12.11) with the mixing matrix 7 ab (iV) (for a,b — q,g) known to two loops for the polarized 
distributions [33] and to three loops for the unpolarized distributions [3]. For our purposes, we 

7 Thc asymptotic series (|2.11l) can be resummcd using the Borel transform 70 (N) = J °° dx e~ x / a j(x), with 
a = l/N being the expansion parameter. Indeed, for z = e~ x this relation coincides with the definition of the 
moments provided that j(x) = Po(e~ x ). 
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have to examine the eigenvalues of the mixing matrix, 7 and 7 9 . They can be found as solutions 
to the characteristic equation 

( la ) 2 - la X(N) + Y(N) = 0, (2.16) 
with the functions X(N) and Y(N) defined as 

X(N) = tr 7 (iV) = 7 W + 7 99 , 

Y (AO = det 7(JV) = 7 99 7 39 - 7 9 V" . (2.17) 

The explicit expressions for roots to (12 . 1 6[) are not particularly suitable for analyzing the large 
iV asymptotics of the anomalous dimensions 7 g and 7 S . As we will see in a moment, it is more 
advantageous to study the properties of X(N) and Y(N). 

As in Sect. 2.1, we first examine the large A" expansion of the one-loop mixing matrix 7q 6 (A"). 
Expanding the well-known expression for ^(N) in powers of the parameter 1/ J defined in (II. 16ft . 
one finds 

7 " = C F [2 In J 2 - 3 + 0(1/ J 2 )] , ^ = -2n f /J + 0(1/ J 2 ), 

(2.18) 

7o OT = -2C F /J + 0(1/J 2 ), 7o" = 2 A" c In J 2 - f3 + 0(1/ J 2 ), 

with (3q given in (12. 5p . The matrix elements ^ are the same for the unpolarized and polarized 
distributions. For 7q 9 their difference scales as ~ 1/J 4 while for 7q 9 and 7q 9 it scales as ~ 1/J 3 . 
The large A" expansion of the diagonal matrix elements, 7^ and 7q 9 , runs in even powers of 1/J 
only. Similar expansion of the off-diagonal matrix elements, 7(j ff and 70 9 , involves all powers of 
1/J, but odd powers of 1/J disappear in their product 7(j ff 7o 9 — ^CpUf/J 2 + 0(1/ J 4 ). This 
implies that the matrix elements ^(N) do not verify the parity preserving relation (11.171) but 
this relation holds true for the trace and determinant of the mixing matrix, X(N) and Y(N), 
respectively. As a consequence, the eigenvalues of the one- loop mixing matrix, Eq. (12.161) . have 
the same property - their asymptotic expansion runs in even powers of 1/J only. E This property 
can be thought of as a generalization of the relation (12.111) for the singlet anomalous dimensions. 

In a close analogy to the nonsinglet case, one can apply fll.l4f) to relate the two eigenvalues 
7 a (N) of the mixing matrix to two different functions f a (N) (with a = q,g). To the lowest order 
in the coupling constant, the functions coincide, 7 a (A r ) = f a (N) + 0(X 2 ), and, therefore, f a (N) 
verify the parity preserving property (I1.17P to one loop accuracy. We will demonstrate in Sect. 3, 
that for the functions f q (N) and f g (N) defined in this way, the property (I1.17P holds true to 
three loops at least. 

2.3. Anomalous dimensions in SYM 

Let us now extend the analysis to twist-two operators in supersymmetric Yang-Mills theories. In 
the simplest case of the SYM theory with N = 1 supercharge, the twist-two operators can be 
obtained from the singlet operators (I2.15P by substituting the quark fields with the gaugino fields 
belonging to the adjoint representation of the SU (N c ) group. As a result, the one-loop anomalous 
dimensions in the M = 1 SYM theory coincide with the QCD expressions for Cp = rif = 
N c [TOl [23]. In the SYM theories with J\f > 1 supercharges, construction of twist-two operators 

8 Provided that the difference 7q 9 — 7q 9 scales as O(J ) which is indeed the case for (|2.18[) since Cp 7^ N c . 



9 



becomes extremely cumbersome due to larger number of fundamental fields involved (gaugino, 
gauge fields and scalars). As was shown in Ref. [M], the problem can be circumvented by 
employing the formulation of SYM theories in terms of light-cone superfields [35]. This approach 
takes a full advantage of superconformal invariance of the SYM theory and provides a unifying 
description of the twist-two (superconformal) operators for different number of supercharges A/. 

Similar to the QCD case, the form of the one-loop mixing matrix in the SYM theory is 
constrained by the (super) conformal invariance of the classical Lagrangian [36]. It allows one to 
classify the twist-two operators according to representations of the superconformal SU(2,2\J\f) 
group and, more precisely, its collinear subgroup SL(2\J\f). In this way, one can show that, in the 
SYM theories with M < 4 supercharges, all twist-two operators belong to two different SL(2\N) 
multiplets whereas in the maximally supersymmetric M = 4 SYM theory they belong to a single 
SL(2\4) multiplet. The reason for such difference is that for M = 4 all 'good' components of the 
fundamental fields, providing the building blocks for the twist-two operators, can be comprised 
into a single light-cone superfield whereas for M < 4 they are described by two different light-cone 
superfields. 

The superconformal symmetry ensures that the twist-two operators inside the supermultiplet 
have the same anomalous dimension. In the SYM theory with M < 4 supercharges, the anomalous 
dimensions of twist-two operators belonging to two supermultiplets are given to one loop by [M] 



7o ,+ (A0 = 4[^(JV + 1)+7e]-A), (2-19) 



7o,- (AO 



.1X^-1)1X4 -AO 



i/)(N + 3 - A") + i;(N - 1) - r (jV + 3 -AT) + ^ 



(3o 



with (3q being the lowest order coefficient of the beta- function defined in f ]2.6j) . In the M = 4 
SYM theory, the twist-two operators belong to the same supermultiplet and their anomalous 
dimension is given to one loop by [37] 



T0l+ (JV)=4[V(iV + l) + 7 E ]. (2.20) 

Here nonnegative integer N parameterizes the superconformal SL(2\M) spin of the twist-two op- 
erators. For the operators with the anomalous dimension 7 0i+ (A) and 7o i _(A r ) the corresponding 
quadratic Casimir is given respectively by 



J2 = N{N + 1) , J 2 _ = (N - fAf) (N + 1 - \M) . (2.21) 

For M = the expressions (12.191) coincide with the one-loop anomalous dimensions of twist- 
two operators in pure gluodynamics. For M = 1 they match the QCD expressions evaluated 
in the supersymmetric limit Cp = rif = N c [101 [23]. In particular, 7 + (A) corresponds to the 
anomalous dimensions of quark transversity and linearly polarized gluon distributions [381 132]. 
For N = even/odd the function 7o ) _(A r ) defines the eigenvalues of the one-loop mixing matrix 
7q & (A) for singlet unpolarized/polarized distributions [23] . 

Let us examine the dependence of the anomalous dimensions (12.191) and (I2.20p on the su- 
perconformal Casimirs (I2.2ip at large N. The expression for the anomalous dimension 7 0j+ (A), 
Eq. (I2.20p . is similar to (12.81) (for rj = — 1) and its asymptotic expansion in powers of can 
be easily obtained from (12.111) . In the SYM theory with M < 4 supercharges, the expression for 
7o,_(A), Eq. (I2.19p . involves the factor (— 1) N and its analytical continuation from even and odd 
N gives rise to two different functions [15] that we shall denote as 70 (A, a) (with a = ±) 

7o (A, a) = 2 [In P + ^( JZ 2 )JZ 2 + <xy? 2 ( J Z 2 ) JZ 4+M ] - A>, (2.22) 
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where, by definition, 7 _(A r ) = jo(N, a = ( — 1)^) for integer N. Here the functions (pi t % are 
given by asymptotic series in 1/J 2 with the M— dependent coefficients. 

We notice that the functions f)2.22|) have different analytical properties for even and odd 
number of supercharges N ' . For M = and M = 2 the functions 70 cr) admit asymptotic 
expansion in integer powers of 1/J 2 only For M = 1 this property is not valid for 70 (iV, cr) 
anymore but it is restored in the sum 70 (iV, +) + 70 (AT, — ). Still, the asymptotic expansion of 
the difference 70 (-/V, +) — 70 (iV, — ) runs in odd powers of 1/ J_ only 

To understand the reason why the parity preserving property is modified for M = 1, we 
recall that jo(N, a) defines the eigenvalues of the singlet mixing matrix in QCD, Eq. ( 12.181) . in 
the super symmetric limit Cp = rif = N c . More precisely, for integer N the solutions to the 
characteristic equation ( 12.161) are given by 71 = 7o,_(iV + 1) and 72 = 7o,-(iV). Then, it follows 
from (I2TT6D and ([2T7D that 

7i - 72 = yj(tf - ll 9 f + 47o"7o 9<? (2-23) 

with the matrix elements of the mixing matrix given by (I2.18p . We demonstrated in Sect. 2.2, that 
(7q 9 — 7q 9 ) 2 and 7o 9 7o 9 are given by asymptotic series in 1/J 2 = 1/(N(N + 1)) and, therefore, 
7i — 72 should also admit similar expansion provided that the leading ~ J° term inside the 
square-root in the right-hand side of (I2.23P is different from zero, or equivalently 71 — 72 7^ at 
large N. In the M = 1 SYM theory, the eigenvalues of the mixing matrix become degenerate 
in the large N limit and the above condition is not satisfied. Indeed, one finds from ( 2.181) that 
(7o 9 — 7o 9 ) 2 + 47o 9 7o 9 ~ 1/ J 2 in the supersymmetric limit, Cp = rif = N c , leading to ^~ 

7o (7V, a) - 7o (iV + 1, -a) = J~ l ip(J~ 2 ) (2.24) 

with J 2 = N(N + 1) and (p(J~ 2 ) given by series in 1/J 2 . For M =1, one replaces 70 (iV, cr) by 
its general expression (12.221) (without assuming parity property of the function ^2), takes into 
account the relation J 2 = (iV— ^)(N+^), Eq. (12 .2 1 j) . and deduces from (I2.24p that, in agreement 
with (12.221) . the cr— dependent term inside 70 (N, a) should admit an asymptotic expansion in odd 
powers of 1/ J_ only. 



3. Parity preserving relations at higher loops 

We demonstrated in Sect. 2, that to one-loop accuracy, the large N expansion of the twist- 
two anomalous dimensions in QCD and in its supersymmetric extensions (except the M = 1 
SYM theory) runs in integer powers of the (super) conformal Casimir J 2 and, therefore, possesses 
the parity preserving property (I1.17p . In the z— representation (11.11) . this translates into the 
reciprocity relation (I1.3P for the corresponding one-loop splitting functions P(z). Going over to 
higher loops, one finds that the conformal spin of Wilson operators get renormalized by an amount 
proportional to their anomalous dimension (II. lip . Moreover, in gauge theories with nonvanishing 
beta-function the anomalous dimensions receive conformal symmetry breaking contributions. We 
will argue in this section, that once both effects are taken into account, the parity preserving 
relation (II. 3p does not hold true for the anomalous dimension 'y(N) but it remains valid (to three 
loops at least) for the function f(N) defined in (I1.14p . 

9 The phenomenon is well-known in quantum mechanics for the system of two almost degenerate energy levels 
7p 9 and 7q 9 whose interaction energy is of the same order as the level splitting, 7q 9 7q 9 ~ (7q 9 — 7o 9 ) ~ l/^ 2 - 
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3.1. Space- like anomalous dimensions 

Let us examine the relation between the space-like anomalous dimension of twist L and the 
function f(N) 

1S (N) = f(N + ^s(N) - \L(3{Xj) . (3.1) 

We recall that the shift of the argument in the right-hand side of (13.11) takes into account renor- 
malization of the conformal spin of the quasipartonic operator of twist L due to nonzero value 
of the anomalous dimension and of the beta- function. The function f(N) admits a perturbative 
expansion in the coupling constant A similar to (12.31) 

f(N) = Xf (N) + X 2 f 1 (N) + ... . (3.2) 

To the lowest order one has fo(N) = 70 (iV). Substituting (13.21) into (13.11) and matching the 
coefficients in front of powers of A, one can establish the relations between higher order corrections 
to f(N) and r ys(N). It is convenient to introduce two functions 

ls(N) = 1S {N) - \L(3{\) , 

f(N) = f{N)-\Lp{\), (3.3) 
so that the relation (13. ip simplifies as 

ls(N)=f(N+\7 S {N)) . (3.4) 
Solving this functional relation one obtains its formal solution as (Lagrange-Burmann formula) 

WW = E h (i^)*" 1 [f(N)] k = f(N) + \ (?{N))' + ^(fW)" + 0{f\N)) , (3.5) 

k=l 

where prime denotes a derivative with respect to N. To determine the function f(N) one has to 
invert this relation. To this end, one uses (13.41) to get 

f(N)=y 8 (N-±f(N)) (3.6) 

and obtains solution to this equation as 

00 1 11 
faO = E u ("l^)"" 1 fis(N)] k = ls(N) -jijUN))' + ^(l 3 s(N))" + Otft(N)) . (3.7) 

k=X 

This relation allows one to reconstruct the function f(N) from known expression for the space-like 
anomalous dimensions 75 (N). 

3.1.1. Parity preserving relations 

We are now ready to verify the parity preserving relation (11.171) for the scaling function f(N) 
in higher loops. For the nonsinglet anomalous dimensions, the analysis goes along the same 
lines as in Sect. 2.1 - one replaces 7,s(A r ) in (13.71) and (13.31) by available two- and three-loop 
nonsinglet anomalous dimensions, evaluates the function f(N) to the same loop accuracy and, 
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then, examines its asymptotic expansion in J 2 = N(N + 1) at large N. For the singlet anomalous 
dimensions, one uses two solutions to the characteristic equation (12.161) . 7 g (A r ) and j g (N), and 
defines the corresponding functions f q (N) and f g (N) with a help of (13. 7p (see Appendix for more 
details). 

In our analysis, we used expressions for multi-loop anomalous dimensions of various twist-two 
operators in QCD and in SYM theories available in the literature. They include 

• Two-loop longitudinally polarized singlet distributions in QCD [33] ; 

• Two-loop gluon linearly polarized distribution in QCD [32j ; 

• Two-loop quark transversity distribution in QCD [38] and its analogs in SYM theories with 
N = 0, 1, 2 supercharges [27] ; 

• Three- loop nonsinglet unpolarized distributions in QCD [2]; 

• Three-loop singlet unpolarized distributions in QCD [3]; 

• Three-loop 'universal' distribution in M = 4 SYM [39J. 

These anomalous dimensions are given by complicated expressions involving nested harmonic 
sums and various color factors. Going over through a lengthy calculation, we worked out their 
large N expansion in powers of J -1 = (N(N + 1)) ^ 2 and split the corresponding series as 

7s(A0 = 7+(ln J, J- 1 ) + j-Qn J, J" 1 ) + |L/3(A) . (3.8) 

Here the beta-function term is introduced for the later convenience and j± are series in J~ l of a 
definite parity, 7±(ln J, — J" 1 ) = ±7±(ln J, J -1 ), 

7+ = r (inJ) + r 2 (inJ) j- 2 + o(j- 4 ), 

7_ = ri(lnJ) J^ + TaQnJ) J~ 3 + 0( J" 5 ) , (3.9) 

where r^ln J) are given by series in the coupling constant with the coefficients depending on 
In J = In J + 7 E . Matching (13. 8p and (jl.7p one can establish the correspondence between the two 
sets of coefficient functions (for L = 2) 

T = A(A)lnJ + B{X) 

T 1 = C_(A)ln J + D_(\). (3.10) 

In all cases except the M = 4 SYM, the resulting multi-loop expressions for the coefficient func- 
tions I\.(ln J) are too cumbersome to be displayed here. For the three- loop 'universal' anomalous 
dimension our analysis gives (see Appendix for details) 

r = 4alnJ-6a 2 C3 + « 3 (-|vr 2 C 3 + 20C5) + C(« 4 ), 

Ti = 8 a 2 lnJ-12a 3 C 3 + C(a 4 ), 

T 2 = |a + 2a 2 + a 3 (|tt 2 + (16- fvr 2 ) lnJ-8 (lnJ) 2 ) +C(a 4 ), 

T 3 = a 2 (|-|lnJ)+« 3 (4+|C3-8 1nJ)+0(« 4 ), ... (3.11) 
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where a stands for the so-called physical coupling constant [40] 

a = ir cusp (A) = A - | vr 2 A 2 + ^ vr 4 A 3 + 0(A 4 ) . (3.12) 

Replacing 7s(iV) in the relation (13.71) by its expression, Eqs. (I3.8p . (13.91) and (13.111) . one finds the 
scaling function in the M = 4 SYM theory after some algebra as 

/(V=4) (JV) = 4aln J_ 6a 2 C3 + a 3(_4 7r 2 C3 + 20C5) 

+ [la + 2a 2 +(l7r 2 -ln 2 \nJ)a 3 ]J- 2 + 0(J-\a A ) . (3.13) 

We observe that, in agreement with (I1.17p . this series does not involve odd powers of 1/J while 
they are present in the series for the anomalous dimension (13.81) . 

We repeated similar analysis for the above-mentioned two-loop and three-loop anomalous 
dimensions of twist two and found that in all cases the corresponding functions f(N) verify 
the parity preserving relation (11.171) ! The properties of the functions f(N) are discussed in the 
Appendix. 

3.1.2. MVV relations 

As was mentioned in the Introduction, the parity preserving property (I1.17P implies certain 
relations between different terms in the large N expansion of the anomalous dimensions. To 
obtain these relations, one notices that according to (11.171) the large N expansion of the deriva- 
tives (c?Ar) n / = (y/T+ l/(4J 2 )dj) n f runs in odd/even powers of 1/J depending on parity of n. 
Matching (I3.5P and (13. 8p one finds that the series 7+ and 7_ are given by the sum of terms 
involving derivatives of the even and odd order, respectively, 

oo 

7+ = Y:^^(ld N ) 2k [f(N)] 2 ^ = f(N) + 

oo 1 

7- = E^l^r -1 ^^ ( 3 - 14 ) 

Inverting the first relation and substituting the resulting expression for f(N) into the second 
relation, one can express 7_ in terms of 7 + . To four- loop accuracy one finds 

7- = l(7^) , + ^(-7 + (7^)' / + i(7t) /, ) / + --- (3-15) 

where prime denotes a derivative with respect to iV and ellipses stand for terms involving higher 
number of 7 + and their derivatives with respect to N. It follows from (13.151) that in order to 
determine 7_ to four loops it is sufficient to substitute 7 + by its three-loop expression. 

Replacing r y± in (13.151) by their expressions ( I3.9P and comparing the coefficients in front of 
odd powers of 1/ J one obtains an infinite system of relations between the coefficient functions 
Yk entering ( 13. 9ft 

Ti = |r r , 

r 3 = -^f + ir2(f ) 2 -ir (r ) 3 + ^roro + |(r r 2 ) -r r 2 , (3.16) 
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where dot denotes a derivative with respect to In J. Here, in the second relation, we took 
into account that To, Eq. f 1 3 . X 1) . is linear in In J and therefore Tq = 0. It is straightforward 
to derive from (13.151) similar relations for the subleading functions r 5 ,r 7 , . . . but they become 
more cumbersome. Substitution of (13.101) into the first relation in (13. 16[) yields C_ = \A 2 and 
D_ = ^A(B — (3), in agreement with the exact three-loop result in QCD, Eqs. (11.91) and (jl.lOp . 
In the M = 4 SYM theory, the relations (13.161) can be easily verified with a help of (13.111) . 

We conclude from (13.151) and (I3.16P that the coefficients in front of odd powers of 1/ J in the 
large N expansion of the anomalous dimension, Eqs. (I3.8P and (I3.9p . can be expressed in terms 
of the coefficients accompanying smaller even powers of 1/ J to less number of loops. 

3.2. Time-like anomalous dimensions 

Let us now extend consideration to the time-like anomalous dimensions 7t (A r ). They describe 
the scale dependence of the partonic fragmentation functions p2] and, in distinction with the 
space-like anomalous dimensions discussed above, are not related to local Wilson operators. In 
what follows we shall assume, following [5], that the time-like anomalous dimensions satisfy the 
relation 

lT (N) = f{N- | TT (iV) - \L(3{\)) . (3.17) 

with the function f(N) the same as for the space-like anomalous dimension ( 13. ip . The relation 
(13. 1 71) generalizes to higher twists the relation ( 11. 5ft proposed in P, [5] and incorporates the 
beta-function contribution. Inverting (13. 1 Tl) and combining it with the similar relation for the 
space-like anomalous dimension (13.11) . one obtains 

/(AT) = 7 t {N + ±/(iV) + \m\)) = Is (N - f /(AO + \m\)) • (3-18) 

One can exclude /(A") from (13. 18[) and obtain two equivalent functional relations between the 
space-like and time-like anomalous dimensions 

ls(N) = lT (N + ls (N)), 

7t(A0 = ls{N - lT {N)) . (3.19) 

Notice that the beta-function term dropped out from these relations. The second relation in 
(I3.19P takes the same form as (ll.4p but in distinction with the latter it now holds true for 
arbitrary A". 

Following [5], one can apply the relations (I3.19P to verify the reciprocity relations (II. 2p . It 
is convenient to switch in (11 .2p from the splitting functions to their moments (11. ip and examine 
the relations between the anomalous dimensions 7s , (A r ) and 7r(A r ). From the second relation in 
(I3.19P one gets similarly to (13.61) and (13. 7[) 

oo 1 

= E If C"^)*" 1 t?s(N)} k ■ (3-20) 

k=l 

We conclude from this relation that the difference between the space-like and time-like anomalous 
dimensions is given by 

7r (A0 - 7s(A0 = -\{l 2 s(N))' + ~(7§(A0)" - ^aW)™ + Otf a (N)) , (3.21) 
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where prime denotes a derivative with respect to N. It is different from zero starting from two 
loops. Moreover, replacing 75 (JV) in the right-hand side of (13.21j) by its three-loop expression, 
one can evaluate 7r(iV) — 75 (A r ) to four-loop accuracy and, then, translate it with a help of fll.il) 
into the difference of the splitting functions Ps{z) ~ Pt{z). 

Let us now examine the second reciprocity relation in (jl.2p and evaluate the moments (11.11) 
of the difference of the splitting functions entering both sides of (11.21) 

A(N) = - f dz z N ~ 1 [P T {z) + zP s (l/z)\ . (3.22) 
Jo 

By definition (II. ip . the moments of Pt{ z ) give rise to the time-like anomalous dimension 7r(-/V) 
which can be expressed in its turn in terms of the function f(N) using To calculate 

the moments of Pg(l/z) in (13.221) . we go over to the large N limit and change the integration 

variable z = e~ x ^ with j — N + | = (J 2 + |)^ 2 as in (12.131) . Expanding the a:— integral in 
inverse powers of J, one finds that it is given by the asymptotic series (I3.8P with J -1 — > — J -1 

/ dzz N P s (l/z) = 7+ (lnJ,J- 1 ) - 7 _(lnJ,J^ 1 ) + \L(3{\). (3.23) 
Jo 

Then, one takes into account (13.141) and obtains 

00 

A w = E k^-\ d *) k ~ l { [fw + ^( A )] fc - W N ) - \ L ^ x )] k } ■ ( 3 - 24 ) 

k=2 

It follows from this relation that A(iV) vanishes in a conformal gauge theory with /9(A) = 0. 
Otherwise, A(N) receives corrections starting from two loops. 

The relation ( 13. 24ft takes a simple form when expressed in terms of the time-like anomalous 
dimension 7^(A r ). Inverting ( 13.17P and substituting f(N) in ( 13.241) by its expression in terms of 
7t(-/V), one gets after some algebra 

A(N) = lT (N) - lT (N + lLp(\)) . (3.25) 

Its expansion in powers of the beta function reads 

A(N) = [i T {N) + \m\)i+{N) + ±(Lf3(\)) 2 ^(N)] + 0(\ 5 ) . (3.26) 

We observe that this relation has the same remarkable property as Eqs. (I3.15P and (13.211) - using 
the known result for three- loop anomalous dimensions, one can determine A(N) to four loops. 

Our analysis in this section relied on the parity preserving property ( 11. 17ft of the scaling 
function f(N). We verified this property in QCD and in SYM theories using known results for 
the multi-loop anomalous dimensions. In the next two sections, we will argue that the property 
( I1.17P holds true to all loops in QCD in the so-called large (3q limit and in the M = 4 SYM theory 
for special class scalar operators of higher twist in the strong coupling regime. 

10 One observes that the relation (|3.17p transforms into (|3.1| under substitution 7s (AT) — > —jt(N) and f(N) — > 
—/(TV). The expression for jt{N) can be easily obtained by applying the same transformation to (|3.5|) and (13.31) . 
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4. Parity preserving relations in the large — limit 



The anomalous dimensions in QCD depend on the number of quark flavours rif. One can make 
use of this fact to obtain the all-loop expression for the anomalous dimensions j(N) and the 
functions f(N) in the large flo limit 

a = = fixed , @ — > oo . (4.1) 

Here (3q = \-N c ~\ n f * s the one-loop coefficient of the beta-function in QCD and rif is assumed 
to take large negative values. Notice that this limit does not exist in SYM theories, since flo is 
uniquely fixed there by the number of supercharges, Eq. (12. 6[) . 

Let us first consider the nonsinglet (unpolarized) twist-two anomalous dimension in QCD. To 
one-loop order, it takes the form 7 ns = X'Jq(N) + 0(X 2 ) with jo(N) given by (I2.8P for r\ — 1. At 
higher loops, it receives corrections ~ X(Xrif) n enhanced in the large (3q limit. They come from 
the so-called single renormalon chain Feynman diagrams which can be resummed to all orders in 
the rescaled coupling constant a, Eq. (14.1 1) . Then, the nonsinglet anomalous dimension in QCD 
is given in the large (3$— limit by [4TI |4*2] 



^\N) = 2A{a) 



i[)(N + a) -if>(l + a) + 



N-l I a 2 + 2a-l (1 + a) 2 



(l + a)(N + a) (a + 2) (N + 1 + a) 



(4.2) 

where the superscript (oo) indicates the leading asymptotic behaviour for (3q — ► oo and A(a) is 
the cusp anomalous dimension in the limit ( 14.11) 

2C f sin(7ra)r(4 + 2a) 
A[a) ~ 3/3 vrP(a + 2) " (43) 

To lowest order in a, 7^' ) (A r ) coincides with the exact one-loop result (12. 8p (for rj = 1). 

Let us substitute (14. 2p into the relation (13 .ip and reconstruct the corresponding function 
p°°'(N). One finds that the relation (13.11) simplifies significantly in the large— (3q limit because 
the anomalous dimension scales as ^^(N) ~ 1/A) and, therefore, it can be safely neglected in 
the right-hand side of (13. ip . Together with f3(\) = — /? A + . . . = —2a + . . . this leads to 

l { n T ] (N) = f^(N + a). (4.4) 

As a result, one finds from (14. 2p after some algebra 

(2a + 1 + 7]f 3 - a 2 



f^\N) = 2A{a) 



IN(N + 1) 2 (1 + a) (2 + a) 



(4.5) 



with 7] = 1. The reason why we displayed the rj— dependence in (14. 5p is that for 77 = — 1 
the relations (14. 5 p and (14.41) describe the quark transversity anomalous dimension, Eq. (12. 8p . 
in the large p limit [13]. As expected, the function p°°'(N) also depends on the parameter 
e CT = (3(\)/2 = —a + . . . defining the 'critical' value of the space-time dimension d cr = 4 — 2e CT for 
which the gauge theory becomes conformal. Since 'ym(N) ~ l//?o, one deduces from ( I3.19P that 
the space-like and time-like anomalous dimensions coincide in the large (3q limit and, therefore, 
the first of the two reciprocity relations in (11.21) is exact. 
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We recall that the expression (14.511 resums perturbative corrections to the function f m (N) of 
the form ~ A(An/) n . The N— dependence of the function /nf^(iV), Eq. ( 14.5j) . is similar to that of 
the one-loop expression (12.81) - higher order corrections in a only modify the coefficients in front 
of ip(N + 1) and 1/(N(N + 1)). Therefore, in agreement with (11.171) . the large N expansion of 
fnT\N) runs in integer negative powers of J 2 = N(N + 1) only. 

Let us now examine the large (3q limit of the twist-two (un)polarized singlet anomalous dimen- 
sions, 7 g and j g , defined as solutions to the characteristic equation (12.161) . To one-loop accuracy, 
the mixing matrix is given by (12.181) and the nj— dependence resides in two matrix elements only, 
7o 9 and 7q 9 . Then, it is easy to see that eigenvalues of the one- loop mixing matrix scale in the 
large (3q limit ( 14. ip as 

^°°\N) ~ A = a/A) , ^(N) = -AA) = -2a . (4.6) 

Going over to higher loops one finds [HJ Hlj that leading corrections to both anomalous dimen- 
sions take the form ~ A(An/) n ~ a n+1 / ' (3q. The scaling behaviour of the smallest eigenvalue 
7q°°' ) (iV) is the same as of the nonsinglet anomalous dimension 7ns°^(iV). Therefore, in the large 
(3q limit the relation (13.11) between r ) q (N) and the corresponding function f q (N) takes the same 
form as (14.41) 

7;°° ) (iV) = / ( ; oo) (A r + a). (4.7) 

In distinction with 7g°°' ) (iV), the gluon dominated anomalous dimension 7g°°' ) (iV) approaches a 
finite value ( 14.61) and the defining relation (13. ip looks in this case as 

j(°°\N) = /j~> (N + H°°\N) + a) . (4.8) 
Substituting 7g°°' ) (A r ) = —2a + 0(1/ (3 Q ) into this relation, one finds in the large f3 limit 

7f ) (^) = /i 00) W, (4-9) 

with fg°°\N) + 2a = 0(1/ (3q)- The relations (14 .7p and (14. 9p hold true for the unpolarized and 
polarized anomalous dimensions although the corresponding scaling functions f q (N) and f g (N) 
are different in two cases. 

The resummed (un)polarized singlet anomalous dimensions 7g°°' > (iV") and 7g°°' > (iV) have been 
calculated in Refs. [4T| |4~4"]. Their expressions are lengthy and to save space we do not display 
them here. Making use of (14.71) and (14. 9p . it is straightforward to calculate the corresponding 
functions f q ,(N) and fg(N) and, then, work out their large N expansion. In this way, we 
found that for the polarized and unpolarized anomalous dimensions the expansion of f q ,(N) 
and ft\N) runs in integer negative powers of J 2 = N(N + 1) only! We conclude that the 
parity preserving property (I1.17P holds true in QCD for the twist-two (non) singlet anomalous 
dimensions to all loops in the large {3 limit. 

5. Parity preserving relations in the AdS/CFT 

In this section, we will employ the AdS/CFT correspondence to determine the function f(N) in 
the maximally supersymmetric M = 4 theory in the strong coupling regime. 

11 Notice that the expression for the quark dominated polarized anomalous dimensions, Eq.(5.3) in the second 
reference in [5T], contains a misprint. The factor (2/j, + n — 1) should be replaced there by (2/j, + n — 2). We would 
like to thank J. Gracey for clarifying this point. 
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The gauge/string duality allows one to establish the correspondence between the anomalous 
dimensions of Wilson operators in the planar M = 4 SYM theory at strong coupling and energies 
of strings on the AdSsxS 5 background [TJ H5J ED]- This relationship takes a very precise form 
for Wilson operators built from a complex scalar field Z and the light-cone component of the 
covariant derivatives D + = (n ■ D) (see review jl6] and references therein) 



O { L k} (0) = tr D^Z(0)D^Z(0) . . . D k + L Z{0) 



(5.1) 



with ki being nonnegative integer. In gauge theory, O^(0) can be identified as quasipartonic 
operators with twist L and the Lorentz spin iV = £\ k^. The operators (15. ip with the same 
quantum numbers L and N mix under renormalization and their anomalous dimensions can be 
found by diagonalizing the corresponding mixing matrix. For L = 2, the anomalous dimension 
of the operators ( 15. ip is uniquely specified by the total Lorentz spin N whereas for higher twists 
L > 3 the anomalous dimensions occupy a band whose size depends on N (see review |4T] 
and references therein). The AdS/CFT correspondence provides a prediction for the minimal 
anomalous dimension 7z,(iV) of the Wilson operators (15. ip carrying large quantum numbers (twist 
L and/or Lorentz spin N) in the planar M = 4 SYM theory at strong coupling. Namely, 7z,(iV) is 
related to the energy of a (semi) classical folded string which rotates with the angular momentum 
N on the AdS3 part of the target space and whose center-of-mass is also moving with the angular 
momentum L along a big circle of S 5 [301 EI] 

(5.2) 



<y L (N) = E — N — L . 
Here the string energy E is given by the series in 1/ \/~X 



E = y/XEr 



N L \ 



N L \ 



(5.3) 



with Eq being the classical string energy and E\ being one-loop quantum string correction. 

For the folded rotating string, the classical energy E can be found explicitly in terms of the 
Jacobi elliptic IK— and E— functions [31], HH1 H9] . The quantum correction E\ has been recently 
calculated in Ref. [SD] and its expression is more involved. Keeping in ( I5.3P the first term only, 
one finds from (I5.2p the leading order expression for the anomalous dimension as A ^> 1 and 
N/y/X, L/V\ = fixed [SH [52] 



1l(N) = L 



K(t) 6 str -E(r)-^-l 

Ostr 



N/L=- 



E(t) I a str + — 



J str 



(r) ( b stT + — 

Ostr 



(5.4) 



where A' = g 2 N c /(irL) 2 is the so-called BMN coupling [15]. The auxiliary variables r, a str and 
b str parameterize classical rotating string solution and verify the relations 



Ostr = b stj: / yl T , 



J str 



K(r) 



A' 



v_ 

b 2 

"str 



1-1/2 



(5.5) 



Their substitution into ( 15. 4p yields the anomalous dimension ^l(N)/L as a function of N/L 
through parametric dependence of both functions on the auxiliary modular parameter r. 
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Let us apply ( 15. 41) to reconstruct the corresponding function Jl(N). In the M = 4 SYM 
theory, the defining relation (11.121) takes the form 

lL (N) = f L (N + ^ L (N)) . (5.6) 

We recall that the shift of the argument in the right-hand side of this relation takes into account 
renormalization of the conformal spin of the operators ( 15. ip in higher loops. The operators ( 15. II) 
have Lorentz spin N and canonical dimension N + L, so that their 'bare' conformal spin equals 
3 = N+\L. 

We have seen in the previous sections, that for the twist-two operators with large Lorentz spin 

N the function Jl=2{N) admits asymptotic expansion in inverse powers of the SL(2) Casimir 
J 2 = j(j _ !) In t h e 

case under consideration, an important difference as compared with the 
twist two is that, for A 3> 1 and N/y/X, L/\f\ = fixed, both the Lorentz spin N and the 
twist L are necessarily large and, therefore, discussing large N asymptotics of one has to 

distinguish two different limits N ^> L ^> 1 and L > > 1. We will see in a moment that 
the function fi{N) has the parity preserving property in the former limit only. Introducing the 
scaling variable a = N/L (not to be confused with the physical coupling (I3.12p ). one finds the 
quadratic Casimir for the operators ( 15. lft as 



J 2 = L 2 (a + \)' 



L 2 {a + \) 2 \l + 0(l/y/X) 



(5.7) 



L{a + \)_ 

We expect that in the scaling limit, A ^> 1 and a = fixed, the function /z,(iV) should admit 
asymptotic expansion in integer negative powers of J 2 , or equivalently in even negative powers 
of (a + |). To verify this property, we shall apply Eqs. (15 .4p - (15.61) to reconstruct the large N 
expansion of f(N) and, then, examine the resulting series for different values of a = N/L. 

The coupling constant enters into the right-hand side of both relations in (15.41) through the 
BMN coupling A' = q 2 Nr/(irL) 2 . Assuming that the modular parameter r admits a regular 
expansion in powers of A' o 

t = t (°) + A't (1) + ... , (5.8) 
one finds from the second relation in (I5.4p for A' = 

a + - = v ' . (5.9) 

2 2 y/i - r (o)K( r (o)) v ; 

Examining subleading ~ (A') fc corrections to (15. 4p one can express coefficients in terms of 
the leading one r^°\ Then, substitution of (15.81) into the first relation in (15.41) yields asymptotic 
expansion of the anomalous dimension in powers of the BMN coupling pi 1ST] 3 



Jl(N) = L 
The first few coefficients are given by 



A' T (o) + (y) 2 7 (i) + 



(5.10) 



1 



(0) 



7 



(i) 



2 1 
1, 



T 



[(2 



T 



(0) 



)K(t 



(0)> 



2E(r (0) 



4(2 



-(0) 



-(0) 



(r 



(0)\2 



gy/l _ r (o)E(r 



(0)^ 



(5.11) 



12 As was shown in Refs. [5TJE2], this assumption is only justified for X'hi 2 (N/L) < 1. 

13 The one-loop quantum string corrections to the energy of folded rotating string (|5.3p break the BMN scaling 



of the anomalous dimension 50i 
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The relations ( 15. 91) - ( 15. lip define parametric dependence of the anomalous dimension 7i(A r ) on 
the rescaled Lorentz spin a = N/L. The leading ~ A' term in (I5.10p coincides with the one- loop 
expression for the minimal anomalous dimension of the operators (15. ip in the thermodynamical 
limit L > 1 [Ml El]. 

Let us combine together the relations (15.61) and (I5.10p and determine the function fi(N). We 
find that fi{N) also admits an expansion in powers of the BMN coupling 



f L (N) = L [A7 (0) + (A')V (1) + 



(5.12) 



with the coefficients related to (15. lip as 



f(0) 



7 



(0) 



7 



(i) 



;1 



(0) 



WV_I K 4 (r (o ))(r (o)y 



d T ( ) a 



(5.13) 



We observe a remarkable simplification of as compared to 7*- 1 -*, Eq. (15.111) . Let us examine 
the expressions for 7z,(-/V) and for different values of a = N/L. 

L>iV> 1 

In this limit, for a < 1, one finds from (EE} that a = ^(r (0) ) 2 + ^(^Y + ■ • • with r(0) ^ L 
Expanding (I5.10p and (15.121) in powers of r^ ^ one gets the leading term 

1 



7 (o) = f{0) 

and the first subleading correction 

T (!) = - 

;(D = _ 



7T 



-- 7T 
8 

1 ., 



1 2 3 3 

a — -a + -a + . . . 

2 8 



11 



(5.14) 



4 



7T 



a + 2 a 2 a 3 

4 

1 9 3 n 

a + -a - -a 4 
4 8 



(5.15) 



We observe that the two functions, 7x(-/V") and /i(iV), vanish for a — > 0. Expansion of 
involves all powers of a = A^/ L and does not reveal any symmetry!^*! 



iV> L> 1 

In this limit, for a ^> 1, one finds from ( 15. 9p that — > 1. Denoting r 
from (15. 9p for z — >• 

1 2 f 2 — 2 £ + 1 



a + 



2£ 2 



+ O (z 4 ) 



(°) = 1 — 16z 2 one obtains 



(5.16) 



with £ = ha(l/z). One notices that the series inside the square brackets in the right-hand side of 
(I5.16P only involves even powers of z. Therefore, inverting (15.161) one obtains z as a series in odd 
negative powers of (a + |) with the expansion coefficients depending on In (a + |). 



14 Thc situation is similar to that for the one-loop anomalous dimension of twist-two operators in the Af = 4 
SYM theory, Eq. (|2.20p . The anomalous dimension 7o(A^) = 4[ip(N + 1) — ^(1)] vanishes for N = 0, its expansion 
around N = involves all powers of N but its large N expansion runs in integer powers of l/(N(N +1)). 
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Replacing = 1 — 16 z 2 in (15. lip and (15. 13ft and expanding the resulting expressions around 
z = one gets the leading term 



7 



(0) 



/ 



(0) 



-2 
~2T 



+ 4- 



-£ + 1 



+ 0(z*) 



as well as the subleading corrections to the anomalous dimension 'Jl(N), Eq. ( I5.10p . 

1 2( 



7 



(i) 



7 



(2) 



' 1 £ + 1 (3 £ - 5) 
16 £-1 ~ 2{£ - 1) ' 

and to the function / L (iV), Eq. (151211 . 



4£ 4 



+ O {z 3 ) 
43£ 3 + 137£ 2 



(5.17) 



(5.18) 



160^ + 61 



2£ (£ - ly 



+ 0(z 3 ) 



1 , 2(1 + 1) ^ 



+ 1 £ 3 + £ 2 



- 3 



(5.19) 



16 I 



21 {£ - 1)' 



z 2 + O (z 4 ) 



Here as compared with (15. lip and (I5.13P we also included the next-to-subleading correction to 
the functions (QUI) and ([5712]) . 

As follows from (15.171) . expansion of the leading order term, 7^ = runs in even powers of 
z ~ (a + |) _1 . Going over to subleading terms, Eqs. ( 15. 18ft and (I5.19p . one finds that 7^ and 
(with fc = 1,2,...) coincide at the level of ~ z° corrections but deviate otherwise. Remarkably 
enough, the functions do not contain odd powers of z while they are present in 7^. Together 
with (I5.16P this implies that, in agreement with our expectations, the large a expansion of the 
function /l(AT) runs in even negative powers of (a + |) only. 

The relations (I5.18P and (15.191) were derived under tacit assumption that the coefficient func- 
tions entering (I5.10p and (15.121) are uniformly bounded functions of a = N/L. Examining the 
expressions ( I5.17p . (15.181) and (I5.19P one finds that for N ^> L, or equivalently £ — > 00, this 
assumption is justified provided that X'£ 2 1, or equivalently (X/ L 2 )\n 2 (N/ L) I. In the 
opposite limit, for X'£ 2 ^> 1, one finds that the leading corrections ~ (X'£ 2 ) n to (15.121) can be 
resummed to all orders [51] 



Vi + X'i 2 - 1 



g 2 N c 



7i 



lnN + 



(5.20) 



where A' = g 2 N c /(iiL) 2 and ellipses denote terms suppressed by powers of £ ~ In (a + |). We 
notice that despite the fact that each term of the expansion is suppressed by inverse power of L, 
the resummed expression does not depend on the twist L and exhibits a well-known logarithmic 
scaling [30]. For the twist-two anomalous dimensions, this scaling behaviour arises for iV 3> 1. 
Specific feature of operators of higher twist L ^> 1 is that the logarithmic scaling (I5.20p sets up 
for much larger values of the Lorentz spin such that 



A 



N 



— In" — > 1 



L 2 



L 



(5.21) 



In this kinematical region, analysis becomes more involved due to a complicated form of the 
defining relations (15 ,4p and (15.51) . 
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There is however a much simpler way to understand analytical properties of the function 
Jl{N). Let us consider the relation ( 15. 6ft and rewrite it using (15.41) as 



ft 



x{t)L 
2 



) 



L K(r)6 9tr (r)-E(r) 



festr(-r) 



- 1 



(5.22) 



where x is a function of the modular parameter r 



z(t) = E(r) a str (r) - K(r) 



A' 



-1, 



(5.23) 



a s tr(V) 



with a str = b stv /y/\ — r and 6 str defined in (15. 5p . We expect that for large Lorentz spin N 
the asymptotic expansion of the function should run in integer negative powers of J 2 = 

j 2 [l + 0(l/\/\)~\, Eq. (15.71) . with j = N + \L being the conformal spin. Substituting N with 



To verify the parity preserving property (11.171) for fi it is sufficient to check that the right-hand 
side of (15.221) stays invariant under transformation j — > —j modulo '^'-contribution from In j 
terms. Let us show that this transformation corresponds to going around r = 1 in the complex 
r— plane in the right-hand side of (15.241) . According to (I5.16p . the limit N ^> L corresponds to 
t — > 1. Let us introduce a new variable r = 1 — 16z 2 and examine the relations (I5.22p and (15.241) 
for z — > 0. One first verifies that the elliptic functions K(r) and E(r) have logarithmic branch cuts 
that start at r = 1. Their expansion around z = runs in even powers of z 2 with the coefficients 
linear in In z. Examining the relations (15.51) one finds that b s t r has similar analytical properties 
around r = 1, while a stT has additional square-root singularity a str = b str / \/\ — t = 6 str /(4^). 
Therefore, going around r = 1 in the complex r— plane, or equivalently replacing z — > — z, 
one finds that, modulo contribution from logarithmic cuts, the functions K(r), E(r) and 6 str 
stay invariant while a str changes a sign. Going back to (15.221) and (15.241) . we conclude that the 
function fb{N) is indeed invariant under j — > — j and, therefore, it has the parity preserving 
property. 

6. Conclusions 

Anomalous dimensions of Wilson operators with large Lorentz spin N admit an asymptotic 
expansion in inverse powers of N . The leading term of the expansion scales logarithmically 
with the spin and is related to the universal cusp anomalous dimension. The subleading, power 
suppressed, corrections to the anomalous dimensions are not universal and depend on the operator 
under consideration. Motivated by findings of Refs. P, [2} [3j [5], we argued in this paper that the 
subleading corrections satisfy infinite number of coupled relations. They allow one to reconstruct 
the corrections to the anomalous dimensions suppressed by odd powers of J = ^/N(N + 1) 
in terms of the coefficients accompanying smaller even powers of J to less number of loops. 
Remarkably enough, these relations take the same form for different twist-two operators even 
though they relate to each other quantities that are operator dependent. 

We demonstrated that the above-mentioned properties of the subleading corrections natu- 
rally follow from the parity preserving property of the scaling function f(N). This function is 




(5.24) 
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related to the anomalous dimension 7(iV) through the functional relation ( 11.141) which general- 
izes similar relation proposed in P, [5] . It incorporates the constraints imposed by the conformal 
invariance and takes into account the conformal symmetry breaking corrections to the anomalous 
dimensions due to nonzero beta-function. We would like to stress that the functional relation 
fll,14p taken alone does not tell us much about the properties of the anomalous dimensions unless 
it is supplemented by additional condition for the function f{N). The latter is just the parity 
preserving property of the scaling function f(N), Eq. (11.171) . 

It is important to keep in mind that the anomalous dimensions 'y(N) and the scaling functions 
f(N) have (Regge) singularities at negative N, and therefore they cannot be invariant under the 
transformation N — > —1—N, or equivalently J — > — J. The parity preserving property (11.171) only 
holds true for each individual term in the asymptotic series for f(N) and the Regge singularities 
manifest themselves through a factorial growth of the expansion coefficients at higher orders 
in 1/J. It is interesting to note that the parity respecting relation is not a unique feature 
of gauge theories. Using the four-loop result for the twist-two anomalous dimension in scalar 
<fr A — theory [53], we verified that the corresponding scaling function f(N) has the same form as 
(I1.17P with the only difference that the leading In J term is missing. 

In a gauge theory with nonzero beta-function, both the anomalous dimensions and the scal- 
ing functions are scheme-dependent. In our analysis, we employed the dimensional regulariza- 
tion/reduction scheme because the beta- function contribution to the anomalous dimension can 
be described in this scheme by going over to the 'critical' number of space-time dimensions 
d cv = 4 — 2e cr in which the underlying gauge theory becomes conformal [29] . The change of the 
scheme corresponds to a finite renormalization of the Wilson operator, say O^f (0) = Z(N)On{0). 
Then, the anomalous dimension in a new scheme reads 7g R> (./V) = 7s (iV) — /3(\)d\ n \ In Z(N), so 
that the properties of 7s R) (-/V) at large N are now governed by the scheme- dependent factor 
Z(N). Substituting 7^ R) (N) into (13.71) one can reconstruct the corresponding scaling function 
= jf (JV) + 0(P(X)) and verify that it does not satisfy fTTTTTD for generic Z(N). However, 
imposing (11.171) as a condition for Z(N), one can define the renormalization scheme in which the 
subleading corrections to the anomalous dimensions still verify the relations discussed above. 

In this paper, we used explicit expressions for the two-loop and three-loop anomalous dimen- 
sions of twist two in QCD and in SYM theories to verify the parity preserving relation fll . lTj) . 
One may wonder whether the same property holds true for the anomalous dimensions of higher 
twist. Going over to higher twists, one finds that, in general, there are few operators carrying the 
same conformal spin. In distinction with the twist two, the size of their mixing matrix grows with 
N and its eigenvalues can not be written in explicit form as functions of N. This makes the large 
N analysis of higher-twist anomalous dimensions extremely difficult. Tremendous simplification 
occurs for the special subclass of quasipartonic operators due to hidden integrability symme- 
try of their mixing matrix in the planar ('t Hooft) limit. In QCD and in SYM theories with 
M = 0, 1, 2 supercharges, integrability arises for the so-called maximally helicity operators [5^|[54] 
while in the M = 4 theory it gets extended to a larger class of operators [57]. Making use of 
the integrability, one can work out a systematic large N expansion of higher-twist anomalous 
dimensions in the planar approximation. In this way, one finds [581 El] that for the so-called 
single-trace operators built from L fundamental fields belonging to the adjoint representation of 
the SU(N C ) group and carrying the conformal spin s (s = | for scalars, s = 1 for gaugino and 
s = | for gauge strength tensor), the asymptotic expansion of the one-loop anomalous dimen- 

15 Intcgrability phenomenon in four-dimensional Yang-Mills theories has been first discovered in Refs. [541 155] 
in context of high-energy (Regge) asymptotics of the scattering amplitudes. 
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sion runs in integer negative powers of the conformal Casimir J 2 = (N + Ls)(N + Ls — 1). It 
would be interesting to extend analysis to higher orders in the coupling constant and to verify 
the parity respecting relation for the scaling function /r,(iV) in higher loops. The gauge/string 
correspondence suggests that it should hold true at least for the minimal anomalous dimension 
in the strong coupling regime in the M = 4 SYM theory for higher twists and large spins such 
that N > L > 1. 

Note added: When this paper was in writing we learned from Yuri Dokshitzer and Pino 
Marchesini that they studied the reciprocity respecting equation (11.51) in the Af = 4 SYM theory 
and analyzed the properties of the corresponding three-loop scaling function f(N) in important 
kinematical limits. Our results agree with theirs [59] as far as the parity respecting property 
(I1.17P and the structure of the large N expansion of f(N), Eq. (I3.13p . are concerned. 
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A Large TV expansion of the anomalous dimensions 

In this appendix we present expressions for the large N expansion of various twist-two anomalous 
dimensions 7s(iV) and the corresponding functions f(N). They were obtained in collaboration 
with S. Moch. The anomalous dimensions are given by a series (12 .3D in the coupling constant A 
defined in (12.21) . At large their leading asymptotic behaviour is 

ls (N) = 2r cusp (A) In N + O(N ) , (A.l) 

with the cusp anomalous dimension r cusp (A) depending on the gauge theory under consideration. 
Perturbative expansion of 7s (A r ) can be simplified by using the 'physical' coupling [40], a = 
|Tcusp(A). Then, the relations (I2.3P and (13.21) can be rewritten in terms of this coupling as 

7s (iV) = a (4 In N + 7o, P h) + a 2 7i, P h + a 3 72, P h + C(a 4 ) , 
/(TV) = a(21nJ 2 + / ,ph) +« 2 /i,ph + « 3 /2, P h + 0(« 4 ), (A.2) 

where N = Ne^, J = Je 7E and J 2 = N(N + 1) is the conformal Casimir. Here, the 
7ph — coefficients are given by series in integer negative powers of while the / p h— coefficients 
admit expansion in even negative powers of J only, Eq. (I1.17p . 

The multi-loop twist-two anomalous dimensions 7s(A r ) are expressed in terms of generalized 
harmonic (Euler-Zagier) sums. In QCD, they also depend on the SU(N C ) color and SU(rif) 
flavour factors and their explicit expressions are cumbersome [21 E]- In the M = 4 SYM, to three- 
loop accuracy, the perturbative coefficients in (1A.2|) are N c — independent and, as a consequence, 
they are much simpler as compared to the QCD case. For the sake of simplicity, we shall present 
below the explicit expressions f)A.2[) in the M = 4 SYM theory and outline a difference as 
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compared with the QCD expressions. In M = 4 SYM the anomalous dimensions of all twist- 
two operators are given by the same, universal function of A" [39]. To three-loop accuracy, the 
cusp anomalous dimension, or equivalently the physical coupling constant is given by (I3.12|) . 
Expanding the three-loop result of Ref. [39] at large N and matching it with the first relation in 
flAT2|) one finds 

70 , ph = 2N~ l -lN- 2 + ±N- A + 0(N- 6 ) , 

71, P h = -6Cs+ (8 In N)^ 1 + (6 -4 In N)N~ 2 

+ (| lniV-f)iV- 3 + 4iV- 4 + e>(iV- 5 ) , 

72, P h = (-f7r 2 C3 + 20C 5 )-12C 3 iV- 1 

+ (-8 (In N) 2 + (-§ 7T 2 + 16) In N + § tt 2 + 6 Ca) N~ 2 

+ (8 (\nN) 2 + (f tt 2 -32)lnA^+ 12-tt 2 -2C 3 )iV- 3 

+ (-4(lniV) 2 + ^ IniV-f + i|vr 2 )iV- 4 + 0(iV- 5 ) . (A.3) 

One evaluates the corresponding function f(N) with a help of Eqs. (13. 7\i and (13.31) and, then, 
expands it in inverse powers of J 2 = N(N + 1) to get from (1A.2I) 

/o, P h = \3 2 ~ ll ^ 4 + ^ 6 ) > 

A, P h = -6Cs + 2J- 2 + J~ 4 + 0(J- 6 ) , 

/ 2 , ph = (20C5-tvr 2 C 3 ) + (-lvr 2 lnJ+|vr 2 )j- 2 

+ ((6 + |vr 2 )lnJ-|7r 2 -2)j- 4 + 0(j^) . (A.4) 

Examining these relations, one observes that for the anomalous dimension 75 (AQ the coefficients 
in front of negative powers of in the right-hand side of flA.3|) are enhanced by powers of In 
with the leading term 7„ iP h ~ (la N/N J (for n = 1,2). For the function f(N) one finds that, 
in agreement with the parity preserving relation (11.171) . its expansion runs in negative powers of 
J 2 . Remarkably enough, the one-loop and two-loop coefficients, /o, P h and /i iP h respectively, do 
not involve logarithmically enhanced terms but they reappear starting from three loops /2, p h ~ 
In J/ J 2 ~ In N/N 2 (to be compared with 7 2jph ~ (In N) 2 /N 2 ). 

Repeating similar analysis for the twist-two operators in QCD, we found that, in spite of the 
fact that the large expressions for the functions j(N) and f(N) are much more cumbersome, 
the ~ In A" and ~ In J terms inside the 7 P h— and f p h— coefficients have the same properties 
as in the M = 4 theory for the following anomalous dimensions: nonsinglet unpolarized quark 
operators to three loops [2], quark transversity operator to two loops [3H] and linearly polarized 
gluon operator to two loops [32"] . 

For the singlet (un)polarized twist-two anomalous dimensions in QCD, the analysis becomes 
more cumbersome since one has to solve the characteristic equation (12.161) and, then, use its 
solutions, 7 q(A r ) and j g (N), to define the corresponding functions f q (N) and f g (N) with a help 
of (13.71) and (13.31) . As in (12.171) . it is more advantageous to analyze the properties of the functions 
f q (N) + f g (N) and f q (N)f g (N). Making use of ([32]) it is straightforward to verify that 

f q + f g = X-i(x 2 -2f) / + ^(A >3 -3Xl > )" + 0(A 4 ), 

fj g = Y \l-\X' + \(X') 2 + \XX"-\Y"] +(9(A 5 ), (A.5) 
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where the notation was introduced for the functions 



X(N) = % + % = X-2{3(\), 

Y(N) = %% = Y-XP{\) + im}\ (A.6) 

with X and Y given in (12.171) . In the right-hand side of these relations, the N— dependence of 
all functions (except of /3(A)) is tacitly assumed and prime denotes a derivative with respect 
to N. We verified by explicit calculation that for the three-loop singlet unpolarized anomalous 
dimensions [3] and the two-loop singlet polarized anomalous dimensions [33], the right-hand side 
of both relations in flA.5j) has asymptotic expansion in integer powers of 1/J 2 . This implies that 



the corresponding functions f q and f g (or equivalently f q and f g ) possess the parity preserving 
property (11.171) . The resulting expressions can be described as follows. The off-diagonal elements 
of the mixing matrix scale as 7 OT , 7 93 ~ 1/J (see Eq. (12 . 18[) ) and they affect f g and f g at the 
level of 0(J~ 2 ) corrections only. As a result, for both functions, f q and f g , the corresponding 
fph~ coefficients ( 1A.2I) take the form ( 1A.4j) with the only difference that the coefficients in front 
of J~ 2 , J~ 4 , . . . receive an additional contribution given by (an infinite) series in 1/ln J. 
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